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Abstract A combination of diamagnetic pumping and a nonlocal a-effect of the 
Babcock-Lcighton type in a solar dynamo model helps to reproduce observations 
of solar magnetic activity. The period of the solar cycle can be reproduced 
without reducing magnetic diffusivity in the bulk of the convection zone below 
the standard mixing-length value of 10 13 cm 2 s _1 . The simulated global fields are 
antisymmetric about the equator and the toroidal-to-poloidal field ratio is about 
a thousand. The time-latitude diagrams of magnetic fields in the model without 
meridional flow, however, differ from observations. Only when the meridional 
flow is included and the a-effect profile peaking at mid latitudes is applied, can 
the observational butterfly diagrams be reproduced. 

Keywords: Solar Cycle: models • Magnetohydrodynamics (MHD) • Convection 
Zone ■ Turbulence 

1. Introduction 

The aim of this paper is to draw attention to the possibility of resolving sev- 
eral problems of dynamo theory for solar activity by combining diamagnetic 
turbulent pumping and a nonlocal a-effect in a solar dynamo model. 

In the absence of reliable data on magnetic fields in the deep solar interior, 
modeling the solar dynamo remains a controversial issue (for a recent review, see 
ITobiasl[20 09!). Consensus has developed that the dynamo is driven by two basic 
effects. The f2-effect of nonuniform rotation produces a strong toroidal field from 
a poloidal one and the a-effect of cyclonic motion ( |Parker, 1955[ ) regenerates the 
poloidal field. Already the first models of an afi-dynamo produced oscillatory so- 
lutions resembling the solar cycle ( |Leighton, 19 69; Stee nbeck and Krause, 1969] ) . 
If discrepancies between computed and observed parameters within two orders 
of magnitude are tolerable, these pioneering models are quite satisfactory. A 
closer agreement is difficult and may even seem impossible to obtain in view of 
the basic physics of the dynamo process and decades-long practice of dynamo 
simulations. 
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Already IKohlerl (| 1973)) noticed that turbulent magnetic diffusivity should be 
reduced much below the standard mixing- length value of rj T rj 10 13 cm 2 s _1 in 
order to reproduce the observed period of the solar cycle. The time of diffusive 
decay of a magnetic field can be estimated as Td ~ cPr)' 1 , where d is the depth 
of the convection zone. A cyclic dynamo has to regenerate fields in a shorter 
time in order to overpower the diffusive decay, so that P cyc < Td (-P cyc is the 
period of the magnetic cycle). The observed period can be reproduced with a 
diffusivity value much smaller than 10 13 cm 2 s _1 , this case, however, is difficult 
to justify. Eddy viscosity or thermal diffusivity below 10 13 cm 2 s _1 cannot be 
a correct parameterization for solar convection because with such small eddy 
diffusion the external layers of the Sun are still unstable ( Tuominen et al, 1 994; 



Kitchatinov and Mazur, 2000). The assumption of small magnetic diffusion ne- 



eds to explain why the same turbulent mixing transports momentum and heat 
much more efficiently than a magnetic field. Theories of turbulent transport 
coefficients do not support the assumption of large magnetic Prandtl number. 
Direct numerical simulations of Yousef, Brandenburg, an d Rudiger (2003) also 
give this number of the order of unity. 

Observations show a clear predominance of the equator-antisymmetric (dipo- 
lar) component in the global magnetic field of the Sun ( |Stenflo, 1988[ ). The 
symmetric part is relatively small and does not show an 11-year cycle. The 
critical dynamo numbers for the excitation of the global modes of these two types 
of symmetry are usually very close together. The threshold dynamo number for 
the antisymmetric modes may be relatively small, leading to the preference for 
this type of symmetry, but the situation can usually be changed to the preference 
for the symmetric modes by small variations of parameters in a dynamo model. 
In other words, the equatorial symmetry is unstable to small changes in the 
design of dynamo models. Another problem is related to the ratio of toroidal- 
to poloidal-field amplitudes. The ratio is not smaller than a thousand if the 
poloidal field is estimated by its polar value and the toroidal field - by the field 
strength in sunspots. Solar differential rotation of about 30% can produce in the 
11 years of the solar cycle a toroidal field that is at most 40 times stronger than 
the poloidal filed. Strong radial gradient of rotation in the tachocline does not 
change this estimation because the radial field should be as much weaker there 
compared to the meridional field as the radial shear is larger than the latitudinal 
shear (except for in cases of the relic field penetrating from the radiation zone, 
or the field having a structure with multiple latitudinal belts, both cases present 
new problems). 

Another problem of the so-called catastrophic quenching of the a-effect (Gru- 
zinov and Diamond, 1994; [Brandenburg and Subramanian| I2005P is not directly 
related to observations but is 'internal' for theory. It can be outlined as follows. 
The large-scale fields generated by a-effect dynamos are helical. As the magnetic 
helicity is conserved, small-scale magnetic fields attain helicity equal in amount 
and opposite in sign to that of large-scale fields. Helical small-scale fields produce 
their own magnetic a-effect that counteracts the a-effect of whatever origin so 
that the total a-effect strongly diminishes. 

In a recent paper ( Kitchatinov and Olemskoy, 201 la| , we demonstrated that 
the catastrophic quenching is alleviated in a solar dynamo model that combines 
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a nonlocal a-effect and diamagnetic pumping of large-scale fields by turbulent 
convection. In this new publication it is shown that with this model other 
above-mentioned problems also fade. The diamagnetic pumping, which will be 
discussed in the next section, concentrates the magnetic field at the bottom 
of the convection zone. The turbulent diffusion in this region is small. As a 
result, the solar cycle period can be reproduced. The poloidal field in the near- 
base region is much stronger than on the surface and sufficiently strong for 
the differential rotation to generate kilogauss toroidal fields over a solar cycle. 
The critical dynamo number for dipolar dynamo modes is considerably smaller 
compared to quadrupolar modes, so that the model always results in a global 
field of dipolar parity. However, the model still has difficulties in reproducing 
time-latitude diagrams of solar magnetism. We include the meridional flow and 
vary the latitudinal profile of the a-effect to reproduce observational butterfly 
diagrams. 



2. Diamagnetic Pumping 

Diamagnetic pumping of large-scale fields was predicted theoretically long time 
ago ( [Zeldovich, I 957||Radler, I968[ ) but is scarcely known in the dynamo commu- 
nity and is usually ignored in dynamo models (see, however, Riidiger and Bran- 
denburg, 1995; |Kapyla, Korpi, and Tuominen[ 120061 Guerrero and de Gouveia 
Dal Pino, 2008). 

The nature of diamagnetic pumping is well illustrated by the case of in- 
homogencous 2D turbulence with the fluctuating velocity u' uniform along a 
direction defined, say, by the unit vector e and inhomogeneous in the direction 
perpendicular to e: V(m' 2 ) ^ 0, (e ■ V)u' = 0. For this case, IZeldovichl (II957P 
found that the mean electromotive force, £ = (u' x B') (B' is a fluctuating 
magnetic field), can be written as 

£ = -Vx(r? T B), (1) 

if the large-scale field B is also perpendicular to e; r/ T is the turbulent magnetic 
diffusivity. This means that the large-scale field is transported with the effective 
velocity 

U dia = -V?7 T . (2) 

The diamagnetic pumping expels the field from the regions of relatively high 
turbulence intensity. However, for the case of the mean field parallel to the 
direction e, the mean electromotive force reads 

£ = -r/ T V x B (3) 

and there is no diamagnetic pumping. The magnetic field parallel to the e- 
dircction, e||B, behaves like a scalar field. The diamagnetic pumping, therefore, 
is related to the vectorial nature of magnetic fields, and there is no counterpart 
of diamagnetic pumping for scalar fields. 
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Figure 1. Pictorial explanation of diamagnctic pumping (see text). Turbulent intensity in- 
creases from left to right. Undisturbed magnetic lines are shown by dashed lines. When a flux 
tube is displaced by motion, its right boundary is moving faster on average. Displacements to 
the left compress the flux tube and the field strength in the displaced region increases. Material 
is evacuated from the displaced region as shown by arrows. Displacements to the right result 
on average in a decrease in field strength in the displaced region. 

The pumping effect can be interpreted as follows (Figure [1} . Imagine that 
there is a background field perpendicular to both the direction e and the gradient 
of turbulence intensity. If the turbulent motion displaces a flux tube of such a field 
as the one on the left in Figure [TJ i.e., in the direction of decreasing turbulence 
intensity, the left boundary of the tube is on average moving slower than the right 
boundary. The tube is compressed, and the field strength in the displaced region 
increases. Similarly, displacement to the right decreases the field strength in the 
displaced region. Turbulent mixing produces on average the field transport in 
the direction of decreasing turbulent intensity, i.e., turbulent conducting fluids 
behave diamagnetic. 

In the 3D case, the expression for the effective velocity of diamagnetic trans- 
port changes to ( jKrause and Radler, 1980| ) 

Udia^-^V^. (4) 

The diamagnetic pumping can be very efficient near the bottom of the convection 
zone where the intensity of turbulent convection changes sharply with depth 
( Kitchatinov and Riidiger, 2008 ). 

Doubts have been expressed concerning the efficiency of the pumping in a 
nonlinear regime (Vainshtein and Kitchatinov, 1983:). However, turbulent dia- 
magnetism was recently detected in laboratory experiments with liquid sodium 
( |Spence et al., 2 007). Direct numerical simulations also show downward pump- 
ing of large-scale fields near the base of convectively unstable layers (Tobias et al., 
1998; ITobias et all IDorch and Nordlundl |Ossendrijver et al.\ \WM IMTH IM?2l 
Zigler and Riidiger, 2003). Whether the simulations indeed show the turbulent 
diamagnetism is not perfectly clear. Their results were interpreted in terms of the 
Drobyshcvski and Yuferev (1974) effect of topological pumping (jTobias et al\ 
2001] Dorch and Nordlund, 2001). Anyway, the diamagnctic and topological 



SOLA: sdap.tex; 17 August 2011; 2:07; p. 4 



Solar Dynamo Model 



pumpings are of the same sense (downward) near the base of the convection zone. 
We apply the Equation Q for the effective transport velocity in our dynamo 
model. 

It may be noted that turbulent diffusion and pumping are both reduced by 
rotation and/or a magnetic field. ITobias et al.\ (|200ip and |Ziegler and Rudiger| 
(2003]) observed the rotational quenching effect in numerical simulations. An 
analytical theory of turbulent transport predicts that though the turbulent dif- 
fusion and pumping are modified, the relation between their isotropic parts 
survives under the influence of rotation or a magnetic field ( Kitchatinov, 1988). 



3. The Model 

3.1. Dynamo Equations 

Our dynamo model is based on the mean-field induction equation 

— = Vx(VxB+£). (5) 

The large-scale flow is a superposition of rotation and meridional circulation, 

V = e rsm^/M) + -V x (e^^-) , (6) 
p \ r sin J 

where the usual spherical coordinates are used, VI is the characteristic value of 
angular velocity, / is the normalized frequency of differential rotation, ip is the 
meridional flow stream function, p is density, and e^ is the azimuthal unit vector. 
The mean electromotive force is written as 

£ = -v^Vx yrj-B)+A, (7) 

where the first term on the right stands for both turbulent diffusion and the 
diamagnctic pumping of Equation ((4]) . The second term accounts for the a-effcct 
in its nonlocal formulation ( |Brandenburg, Radler, and Schrinner, 2 008) 

A = y"a(r,r')B(r') dV. (8) 

The local a-cffect is subject to the catastrophic quenching mechanism (Bran- 
denburg and Subramanian, 2005) while the nonlocal a-effect is not (Kitchatinov 
and Olemskoy, 2011a), and for this reason it may be dominating in the dy- 
namo process. Another possibility to avoid the catastrophic quenching is to 
invoke sufficiently efficient fluxes of magnetic helicity (Guerrero, Chatterjee, and 
Brandenburg, 2010). 

Similar to the mean flow of Equation (JB|) , the magnetic field can be written 
as a superposition of its toroidal and poloidal parts, 

B = e^B + V x ( e^— T-— : ) , (9) 
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where A is the poloidal field potential. 

Here we introduce normalized variables. The time is measured in units of 
Rq/tjq; rjo is the characteristic value of the eddy diffusivity. The magnetic field is 
normalized to the field strength Bo for which nonlinear effects become essential, 
and the a-parameter is normalized to its characteristic value ao- The poloidal 
field potential is measured in units of ao-Bo-Ro/^o- The density is normalized to 
its surface value po, and the stream function of meridional flow is measured in 
units of PqR^qVq', Vq is the amplitude of the surface meridional flow. From now 
on, the same notations are kept for the normalized variables as used before for 
their unnormalized counterparts, except for the fractional radius x = r / Rq and 
normalized diffusivity ?y = 7] T /r]Q. The normalized equation for the toroidal field 
reads 

dB r] d ( 1 d{sm6B)\ 1 8 ( d(y/fj xB) 



dt x 2 86 \ sin 9 89 J x dx \ dx 
R m d ( B 8ip\ R m d ( B dip 



x 89 \px sin 9 dx J x dx \px sin 9 d9 
V fdfdA dfdA 



+ 



x \8x 89 d9 dx 



(10) 



where 



V = , (11 



is the dynamo number and 



Rm = ^ (12) 

Vo 

is the magnetic Reynolds number for the meridional flow. The afJ-approximation 
is applied to neglect the alpha-effect in the toroidal field Equation (TTU)) . This 
equation describes the toroidal-field production by differential rotation, its ad- 
vection by the meridional flow, diamagnetic pumping, and turbulent diffusion. 
The poloidal-field equation with nonlocal a-effect is written as 

dA r) . 8 ( 1 dA\ 8 ( ^8A 



dt x 2S ' md M{^^) + ^dx~{ VTl dx~ 
R m ( dtp 8 A dip 8 A 



px 2 sin 9 \8x 89 89 dx 

X 

x sin 9 cos 9 / a(x,x')B(x', 9) dx', (13) 



where x\ is the radius of the inner boundary. The integration in this equation is 
only in the radius with the upper limit x. This qualitatively reflects the fact that 
the nonlocal a-effect at some point x is contributed by the buoyant magnetic 
loops rising from deeper layers {x! < x) and that buoyant velocities are almost 
vertical. Our dynamo Equations (|10|) and (|13[) are very similar to those in our 
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previous publication (Kitchatinov and Olcmsko y, 2011a| , but now we include 
meridional flow and neglect the (local) magnetic a-effect which was shown to be 
insignificant when nonlocal a-effect is allowed for. 

Our boundary conditions assume a perfect conductor beneath the inner bo- 
undary of radius x\ and pseudo- vacuum conditions on the top, 

0, A = for i = Xi, 
0, B = 0, for x = 1. (14) 

The bottom boundary in our model is at x\ = 0.7. 

The initial-value problem for the dynamo Equations (fTO)) and (fl"3|) was solved 
numerically by the grid-point method and explicit time-stepping. The diamag- 
netic pumping leads to a high concentration of the magnetic field near the 
bottom. To resolve fine structures near the bottom, a nonuniform grid was 
applied over the radius with the grid spacing Ax ~ rfl 2 . The grid over the 
latitude was uniform. 

The equatorial symmetry was usually not prescribed. The field was evolved in 
time starting from a mixed-parity initial field and the solution relaxed eventually 
to a certain equatorial symmetry. In order to determine the critical dynamo 
numbers for the excitation of the dipolar (B(0) = —B(tt — 9)) and quadrupolar 
(B(9) = B(ir — 8)) dynamo modes, additional boundary conditions selecting the 
field mode of certain equatorial symmetry were imposed on the equator. 

3.2. Model Design 

3.2.1. Differential Rotation 

For the differential rotation, we use the approximation of hclioscismological data 
suggested by |Belvedere, Kuzanyan, and Sokoloff] (|2000|) 

1 2 4 

f(?> fl ) = ^E cos ( 2m (| - e )) E c ^ n . (15) 

m— n— 

The coefficients C nm of this equation are given in Table 1 of Belvedere, Kuzan- 
yan, and Sokoloff (2000). Figure [5] shows the angular velocity contours. 

3.2.2. Profiles of the a- Effect and Diffusivity 

The kernel function of the nonlocal a-effect in the poloidal field Equation ((15)) 



Xb)/h b )) , 

x a )/h a )), (16) 



9 (y/rjxB) 
dx 

OA 

dx 



was prescribed as follows, 

a(x, x) 



(j>h{x')(j) a (x) 

i + B 2 ( X ',e)' 



cp h (x') - 1(1 -erf ((x 1 
<f> a (x) = -(l + erf((x- 
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where erf is the error function and B 2 in the denominator of the first equation 
accounts for the usual algebraic quenching of the a-effect. We always use x^ = 
X[ + 2.5/ib and x a = 1 — 2.5h a to ensure smoothness of the kernel functions in 
the simulation domain. The parameter /ib-parameter represents the thickness of 
the near-bottom region of toroidal magnetic fields producing the a-effect. The 
parameter h a represents the thickness of the near-surface layer where this effect 
is produced. The a-effect with the kernel function of Equation (fTB]) is very close 



to the nonlocal model of Brandenburg and Kapyla] (|2007|) . It is also similar to 



the Babcock-Lcighton mechanism for the poloidal field production used in the 
dynamo models of |Durne"y| (|1995p and Dikpati and Charbonneau (|1999p . 



The turbulent diffusivity in the bulk of the convection zone varies slightly with 
depth. Near the base of the convection zone, it drops sharply with increasing 
depth. The diffusivity profile of our model reads 

V(x) = r?in + ^(1 - Vin) (l + erf ) , (17) 



where ?7i n is the ratio of diffusivity near the inner boundary to its value in the 
bulk of the convection zone. Computations were performed for the smallest value 
of ?7m = 10~ 4 we were able to apply. 

The major part of our computations was performed with the following values 
of parameters: x n = 0.74, h v = 0.01, h a = 0.02, and hb = 0.002. The profiles 
of diffusivity and kernel functions of the a-effect of Equation ([To]) for this set 
of parameters are shown in Figure [3] Dependence of the results on the model 
parameters hh and h a will be discussed at the beginning of Section @] to explain 
this choice. 
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Figure 3. Profiles of the normalized diffusivity and the kernel functions of Equation (1 1 6 1 of 
the nonlocal a-effect. 



3.2.3. Meridional Flow 

A poleward meridional flow of the order of 10 ms" 1 is observed on the solar 
surface ( jKomm, Howard, and Harvey, 1993] ). Helioseismology confirms that the 



flow persists up to a depth of about 12 Mm ( Zhao and Kosovichev, 2004 ). Theo- 
retical modeling remains the only source of knowledge about the flow in deeper 
regions. Recent simulations predict that one cell of meridional circulation occu- 
pies the entire thickness of the convection zone and the return flow at the bottom 
is not small compared to the surface (Kitchati nov and Olemskoy, 20115] ). Be- 
neath the convection zone, the flow is small (Gilman and Micsch, 2004; Kitchati- 
nov and Riidigcr, 2006). All these findings are qualitatively reflected by the 
following representation for the stream function of the meridional flow 

ip = —cos 9 sin 2 9(j)(x), 

J p{x')^{x'){x' - xs) dx' for x > x s 

(18) 

C j p(x')rf(x'){x s — x') dx' for x < x s . 

Xi 

In this equation, x s is the radius of the stagnation point where the meridional 
velocity changes sign, and C is a parameter whose value is adjusted to ensure 
continuity of the stream function at the stagnation point. 

The stratification of the convection zone is almost adiabatic. An adiabatic 
profile for ideal gas was used for the normalized density, 

p(x)=(l + C p (±-l\ \ ' , (19) 
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Figure 4. Radial profile of meridional velocity in units of Vb/2 at a latitude of 45°. 

where C p = 10 3 . 

The profile of the normalized meridional velocity of our model is shown in 
Figure |U The radius of the stagnation point, x s ~ 0.84, was adjusted so that 
the amplitudes of the flow above and below this point are equal. We set the 
parameter p = 0.25 in Equation (fT8"| . so that the bottom flow is about ten times 
slower than the maximum velocity of the deep equator-ward flow. 



4. Results and Discussion 

4.1. Models without Meridional Flow 

We will discuss the models without meridional flow (R m — 0) first. Then, we 
will include the flow to see what effect it produces. 




10" 4 10" 3 10' 4 10' 3 

Figure 5. Dependence of the critical dynamo number (a) and the magnetic cycle period (b) 
on the thickness h^, of the near-bottom layer producing the cr-effect. Full and dotted lines show 
the results for dipolar and quadrupolar dynamo modes, respectively. 
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Figure [5] shows the dependencies of the critical dynamo number and magnetic 
cycle period on the thickness of the near-bottom layer producing the a-cffcct 
(for constant h t = 0.02). The cycle period P cyc is defined as half the complete 
period of magnetic oscillations in order to compare with the 11-year solar cycle. 
The dipolar parity and the observed period of magnetic oscillations can be re- 
produced only with a sufficiently thin bottom layer, < 4 x 10 -3 . Only in this 
case the equator-antisymmetric modes of the magnetic field are preferred and 
the cycle period is sufficiently long. For the value 770 = 10 13 cm 2 s" 1 of magnetic 
diffusivity, the diffusion time Rq/tjo sa 15 years and the cycle period in physical 
units, P cyc R,Q/r]o, is close to 11 years for small h^. Also for a thin bottom layer 
only, the critical dynamo number for the excitation of magnetic fields of dipolar 
parity is about 30% smaller compared to quadrupolar parity modes so that the 
model produces equator-antisymmetric global fields similar to the Sun. A strong 
decrease of magnetic diffusivity towards the bottom boundary is also impor- 
tant for preference of dipolar parity flChatterjee, Nandy, and Choudhuri7"2 004; 
|Hotta and Yokoyama, 2010 ). We conclude that our model can reproduce a solar 
cycle only if the a-effect is produced by sufficiently deep-seated toroidal fields. 





0.01 0.03 0.05 0.07 
h, 



Figure 6. Dependence of the critical dynamo number (a) and the magnetic cycle period (b) on 
the thickness ht of the near-top layer where the a-effect is produced. Full and dotted lines show 
the results for dipolar and quadrupolar modes of magnetic field, respectively, = 2 X 10 — 3 . 

The dynamo numbers shown in Figure [5] have a minimum at about hy, = 
2 x 10 -3 , and all the results to follow were produced with this value. For smaller 
hh, the critical dynamo numbers are larger because of the decrease in the toroidal 
magnetic flux producing the a-cffcct. For larger h^, the bottom layer is thicker 
than the vertical scale of the toroidal field and the fields of opposite signs partly 
cancel each other's contribution to the a-effect. The dependence on the thickness 
ht of the top layer where the a-effect is produced is simpler. The cycle period 
shown in Figure |5] does not depend on h t until the thickness becomes so large 
that the top and bottom layers overlap. There is also a smooth decrease in 
the critical dynamo numbers with increasing h t . We fix h t = 0.02 for further 
computations. 

Figure [7] shows the magnetic field patterns in the meridional cross-section 
for several instances of a cycle. The magnetic field is highly concentrated at 
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Figure 7. Contours of toroidal-field (top row) and poloidal field lines (bottom row) for several 
instances of a magnetic cycle. The time of the run in units of Rq/t^q is shown at the top. Full 
(dotted) lines show positive (negative) levels and clockwise (anticlockwise) circulation. The 
pictures of the upper row are re-scaled so that the upper (dashed) boundary shows the radius 
of r = 0.74/Jq , below which the toroidal fields arc localized. The dynamo number D = 2.2 X 10 4 
is slightly above the critical value of 1.9 X 10 4 . 



the bottom. The concentration is provided by the diamagnetic pumping. The 
pumping, however, does not affect the (radial) component of the field parallel to 
the pumping direction and the poloidal field lines can come to the surface. The 
toroidal field is confined below the radius of x v where the magnetic diffusivity 
drops sharply with depth. The poloidal field near the bottom is also much 
stronger than at the top. As a result, the toroidal field amplitude in our model 
is about a thousand times stronger than the amplitude of the surface polar field. 
If the poloidal field at the bottom were of the same order as at the top, the 
differential rotation would not be able to produce a strong toroidal field over a 
cycle. Note that the deep toroidal field of several kilogauss in our model is the 
mean field. The total field can include strong fluctuations on the background of 
the mean field. 

Figure [5] shows butterfly diagrams for the surface radial field and deep toroidal 
field. The toroidal field diagram shows the contours of the quantity 



B = sin t 



4>\,{x)B{x) dx, 



(20) 



SOLA: sdap.tex; 17 August 2011; 2:07; p. 12 



Solar Dynamo Model 



90 




-90 I I 

22.0 22.5 23.0 23.5 24.0 24.5 25.0 




22.0 22.5 23.0 23.5 24.0 24.5 25.0 
DIFFUSIVE TIME 



Figure 8. Butterfly diagram of the depth-integrated toroidal field B of Equation (1 20 II (top 
panel) and the surface radial field (bottom) for the model without meridional flow. Time is 
shown in units of R'i./ijo- D = 2.2 X 10 4 . 



to which the surface a-effect of our model is proportional. The factor sin 9 in 
Equation (f2U)) accounts for the dependence of the length of toroidal flux tubes on 
latitude (it is supposed that the probability of sunspot production is proportional 
to the length of the tube). 




Figure 9. Time-latitude diagrams for the model with meridional flow. A standard latitudinal 
profile of the a-effect, a ~ cos 9, was used. The dynamo number D = 2.6 X 10 4 is slightly above 
the critical value of D c = 2.36 X 10 4 . The top and bottom panels show the deep toroidal field 
and the surface radial field respectively. 
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The radial field diagram of Figure [5] is similar to observational diagrams of 
IStenflol ([1988)) and lObridko et al.\ (j2006|) . The toroidal field diagram, however, 
has too broad 'wings' and shows too slow an equatorial drift compared to the 
butterfly diagram of sunspots. There is a possibility to improve the agreement 
with allowance for meridional flow. 

4.2. Models with Meridional Flow 

All computations with meridional flow were performed for the magnetic Reynolds 
number R m = 10 that is a plausible solar value for the magnetic diffusivity 
770 = 10 13 cm s . The 'local' Reynolds number in the near-bottom region of 
low diffusivity is, of course, much larger. It may be expected, therefore, that 
the meridional flow influences primarily the field migration near the bottom. 
Indeed, the toroidal field diagram of Figure shows much faster equatorial drift 
compared to Figure [5] 



90 




-90 I , , , , , I 
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Figure 10. Timc-latitudo diagrams for the model with meridional flow and sin 2 8 cos 0-profilc 
of the a-effect. Computations were performed for the dynamo number D = 4.2 X 10 4 , slightly 
above the critical value of D c = 3.96 X 10 4 . The top and bottom panels show the deep toroidal 
field and the surface radial field respectively. 

Now, however, the poloidal field diagram becomes unsatisfactory. It shows 
equatorial drift also at high latitudes where the poleward migration is actually 
observed ( |Stenflo, 1988] ). The magnetic diffusion near the surface is high and 
the meridional flow does not influence the surface fields directly. However, the 
surface poloidal fields are related to the deep toroidal fields by the a-effect and 
follow the equatorial migration of the deep fields. 

The polar drift of the high-latitude poloidal field can be re-established by a 
change in the latitudinal profile of the a-effect. Observations indicate that the 
current hclicity of solar active regions (the a-effect proxy) docs not increase 
steadily with latitude but has a maximum at mid latitudes (Sokoloff et ai, 



SOLA: sdap.tex; 17 August 2011; 2:07; p. 14 



Solar Dynamo Model 



2008). Theoretical tilt- angles of bipolar magnetic regions show also humps in the 
latitudinal profiles ( |D'Silva and Choudhuri, 1993[ ) . Figure [TU] shows the model 
results for the sin 2 9 cos ^-profile of the a-effect, i.e., sin9 in the last term of 
Equation (TT3"]) was changed to sin 3 9. With this profile, the poloidal field diagram 
is in closer agreement with observations. 

Simulations that take into account the meridional flow also result in magnetic 
fields of dipolar parity. The critical dynamo numbers for the excitation of dipolar 
and quadrupolar dynamo modes in the model with the a-effect peaking at mid 
latitudes (a ~ sin 2 9 cos 9) arc = 3.96 x 10 4 and D°f = 4.64 x 10 4 , respectively. 
In the model where the a-effect peaks at poles (a ~ cos (9), the critical dynamo 
numbers are closer to each other, = 2.36 x 10 4 and D% = 2.49 x 10 4 . Neverthe- 
less, this model also shows the preferred excitation of the equator-antisymmetric 
magnetic fields. 

Though the meridional flow is important for the latitudinal drift of the mag- 
netic fields in our model, the model does not belong to the so-called advection- 
dominated dynamos. Figure [TUJ does not show high concentration of surface 
fields towards the poles typical of advection-dominated dynamo models. This is 
because of relatively large magnetic diffusivity near the surface. The cycle period 
in our model is not controlled by the meridional flow. A sufficiently long cycle 
is realized mainly due to the low magnetic diffusivity in the thin layer near the 
bottom boundary and partly due to the diamagnetic pumping to this layer. 

To probe for the effect of diamagnetic pumping, the computations were re- 
peated with the pumping switched off (for the model with meridional flow 
and mid-latitude peaking a-effect). The cycle period at the critical dynamo 
number reduced to P cyc = 0.47, i.e., a bit less than twice. Other model param- 
eters are influenced much more. Without diamagnetic pumping, the amount 
of toroidal magnetic flux in the near-bottom region, producing the a-effect, 
decreases strongly. Accordingly, the critical dynamo numbers increase by almost 
two orders of magnitude, to = 2.47 x 10 6 and D% = 2.49 x 10 6 with no 
clear preference for dipolar parity. Without the pumping, the poloidal field at all 
depths is of the same order as at the surface. Accordingly, the ratio of amplitudes 
of the toroidal field to the surface polar field decreased to about 43 from 920 in 
the model with diamagnetic pumping. 

5. Summary 

Joint application of diamagnetic pumping and a nonlocal a-effect strongly in- 
fluences the solar-type dynamo model to bring it generally closer to observa- 
tions. The pumping concentrates magnetic fields near the base of the convection 
zone where diffusivity is relatively small. As a result, the solar cycle period 
can be reproduced with the standard mixing-length value of eddy diffusivity, 
r) T rj 10 13 cm 2 s _1 , in the bulk of the convection zone. The near-base con- 
centration of the poloidal field makes it possible for differential rotation to 
wind a strong toroidal field over a solar cycle so that the toroidal field am- 
plitude in the dynamo model is about a thousand times stronger than the 
surface polar field. The model produces equator-antisymmetric global fields. 
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The combination of diamagnetic pumping and nonlocal a-effect also resolves the 
theoretical problem of the so-called catastrophic quenching of the alpha-effect 



(Kitchatinov and Olemskoy, 2011b). 

The model without meridional flow docs not, however, reproduce the butterfly 
diagram of sunspot activity. The time-latitude diagram of the toroidal field 
shows too slow equatorial drift and too broad latitudinal distribution in this case 
(Figure [5]). This can be corrected by allowance for the meridional circulation and 
the change of the latitudinal distribution of the a-effect from a pole-peaked to 
a mid-latitudc-pcaked profile (Figure [TU|) . 

The prescription for the a-effect remains the main uncertainty of the model. 
Operation of the nonlocal a-effect o f the Babcock-Leighton t ype on the Sun 



is plausible from both observational (Dasi-Espuig et al., 20101 and theoretical 



(Caligari, Moreno-Insertis, and Schussler, 1995) evidences. The mathematical 
formulation of this effect remains, however, uncertain. If our model is considered 
as an inverse problem, it predicts that the a-effect is produced by deep-seated 
toroidal fields and it is largest at mid latitudes. 
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